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Dedekind cotangent sums 


Matthias Beck 


Abstract: Let a, ai,..., be positive integers, mi,..., nonnegative integers, and zi, ..., complex numbers. We study 
expressions of the form 

e . 

fc mod a j= 1 ' ' 

Here the sum is taken over all k mod a for which the summand is not singular. These Dedekind cotangent sums generalize 
and unify various arithmetic sums introduced by Dedekind, Rademacher, Apostol, Carlitz, Zagier, Berndt, Meyer, Sczech, 
and Dieter. Generalized Dedekind sums appear in various areas such as analytic and algebraic number theory, topology, 
algebraic and combinatorial geometry, and algorithmic complexity. We prove reciprocity laws, Petersson-Knopp identities, and 
computability statements for the Dedekind cotangent sums. 


1 Introduction 


While studying the transformation properties of 


r](z) 


. _ ^niz/12 n (i - e 2ninz 

n> 1 


under SL 2 (Z), Dedekind, in the 1880’s [De], naturally arrived at the following arithmetic function: Let ((a;)) 
be the sawtooth function defined by 


(Cd) 


{a:} — \ if x Z 
0 if x £ Z . 


(1) 


Here {a;} = x—[x] denotes the fractional part of x. For a, b £ N := {n £ Z : n > 0}, we define the Dedekind 
sum as 



( 2 ) 


The Dedekind sums and their generalizations have since intrigued mathematicians from various areas such 
as analytic [De, Dil, Al] and algebraic number theory [Me, S], topology [HZ, Za, MS], algebraic [P, BV] and 
combinatorial geometry [Mo, DR], and algorithmic complexity [Knu]. 


By means of the discrete Fourier series of the sawtooth function (see, for example, [RG, p. 14]), 


. p —l 


II \ \ l \ - / 7 At 


k= 1 


irk 


0 2nikn/p 


it is not hard to write the Dedekind sum in terms of cotangents: 


1 Ei 7T ka 7 rk 

s(a, b) = — > cot —-— cot — 
46 ' 6 6 

fc=l 
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( 3 ) 


( 4 ) 
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Starting with these two representations (2) and (4) of s(a,b), various generalizations of the Dedekind sum 
have been introduced. This paper constitutes an attempt to unify generalizations of the Dedekincl sum in 
its ‘cotangent representation’ (4). Through the discrete Fourier series (3), this includes most generalizations 
of the ‘sawtooth representation’ (2) of the Dedekind sum. 

Let cot^ 1 ™) denote the m’th derivative of the cotangent function. Our generalized Dedekind sum is introduced 
as follows: 

Definition 1 For ao, ..., a<j € N, mo, ..., rri d 6 Mo := N U {0}, Zq, ... ,Zd £ C, we define the Dedekind 
cotangent sum as 



where the sum is taken over all k mod a o for which the summand is not singular. 

The Dedekind cotangent sums include as special cases various generalized Dedekind sums introduced by 
Rademacher [R], Apostol [Apl], Carlitz [Cl], Zagier [Za], Berndt [Be], Meyer, Sczech [MS], and Dieter [Di2]. 
Section 2 contains the definitions of these previously defined sums and the connections to our new definition. 

The most fundamental and important theorems for any of the generalized Dedekind sums are the reciprocity 
laws: an appropriate sum of generalized Dedekind sums (usually permuting the arguments in a cyclic fashion) 
gives a simple rational expression. The famous reciprocity law for the classical Dedekind sum is as old as 
the sum itself: 

Theorem 1 (Dedekind) If a,b eN are relatively prime then 

, ,, , 1 1 fa 1 b\ 

5 (a,b)+s(b,a)--- + -^- + - + -j . 

The following reciprocity law for the Dedekincl cotangent sums will be proved in section 3: 

Theorem 2 Let ao, •.., a d £ N, mo,... mj £ No, Zq, .. ., z d £ C. If for all distinct i,j € (0,..., d} and all 
m, n £ Z, 

m + Zi n + Zj 

-- , 

CLi dj 

then 

d l Q l l d ( CLn CL 0 * * * CL n ' * * ^d \ 

E l -t\m , a 0 ‘“ a n'“ a d , '- 1 , \ 

(-1) m n \ 2^ —-p-— c I m n mo + lo ••• m n j+l n ■■■ m d + ld I 

"=° (o.-.jin,..., i d >0 <0 ■■■■In-■■■Id.- \ Zn Zq Z n ••• Zd ) 

-- d = rn n 

( (—l) d / 2 if all mfe = 0 and d is even, 

\ 0 otherwise. 


As usual, x n means we omit the term x n . 
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In section 4, we obtain as corollaries the reciprocity laws for Dcdekind’s, Apostol’s, Meyer-Sczech’s, Dieter’s 
and Zagier’s generalized Dedekind sums. The purpose here is not so much a rederivation of old theorems 
but rather to show a common thread to all of them. 

Another basic identity on the classical Dedekind sum is about a hundred years younger than Dedekind’s 
reciprocity law [Kno]: 


Theorem 3 (Petersson-Knopp) Suppose a, b £ N are relatively prime. Then 

yy 5 (J—^b + ka, acfj = a(n) s(b, a) . 

d\n k mod d 

Here a(n ) denotes the sum of the positive divisors of n. 


( 5 ) 


This result has also been extended to certain generalized Dedekincl sums [AV, PR, Zh]. The respective 
identity for the cotangent sum follows from a much more general theorem, stated and proved in section 5: 


Theorem 4 For n, ao, ..., ad € N, mo,..., m d € No, 

"y ' ^mo + l—mt - md—d 


a 0 b 

c I m o 

ri,...,rd mod b \ 0 


a 0 

— Ti (7— rni -rad —1 (ji) c ( rno 

0 

Here a m {n) := Ed|n rfm - 

Finally, we show in section 6 a computability result: 


a i 
m i 

0 


input size of ao ,..., ad- 


fax +ria 0 
0 

• a d 

■ m d 

■ 0 


( a 0 

ax ■ ■ 

• a d \ 


mi 

• rn d 1 

\ zo 

Zx • ■ 

z d ) 


fa d + r d a 0 
m d 
0 


Computability for any of the generalized Dedekind sums involving two integer arguments follows easily from 
the two-term reciprocity laws. However, the reciprocity laws for ‘higher-dimensional’ analogues do not imply 
nice computability properties. 


2 Various Dedekind sums 


In this section we will give an overview of previously defined generalizations of the Dcdekind sum. We do 
not claim any completeness but hope to give a good picture of what has been introduced in the past. 

The sawtooth function ((x)) defined in (1) is the first Bernoulli function Bx(x ), defined to be zero at the 
integers. We used the slightly antiquated function ((a;)) partly for historical reasons, partly because it has a 
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simpler discrete Fourier series. The Bernoulli polynomials B k {x) are defined through 


The first few of them are 


ze 


e— 


_ Bk(x) ~ k 

1 ^ Id 

k> 0 


( 6 ) 


Bi(x) = x- - 

B 2 (x) = x 2 - x + \ 

6 

„ , . 3 3 2 1 

B 3 (x) = X 6 - -X 2 + —x 


1 


B±{x) = x 4 — 2x 3 + x 2 — 7 ^- 

oU 

B 5 (x) = x 5 - ^x 4 + ^x 3 - ^x 
B e{x ) = x 6 - 3x 5 + ^x 4 + \ x2 + ^ 
B 7 {x) = x 7 - \x 6 + E 5 + \x 3 + \x . 


6 


6 


The Bernoulli numbers are B k := B k { 0). The Bernoulli functions Bk(x) are the periodized Bernoulli 
polynomials: 

B k (x) := B k ({x}) . 

Apostol [Apl] replaced one of the sawtooth functions in (2) by an arbitrary Bernoulli function: 


k mod b 


£ tMt) ■ 


( 7 ) 


Apostol’s idea was generalized by Carlitz [Cl] and Mikolas [Mi] to what we would like to call the Dedekind 
Bernoulli sum, defined for a, b , c, m, n £ N as 


^ n (a;b,c):= V] B m (—'] B n (—^\ . 

k mod b V a / V a J 


Another way of generalizing (2) is to shift the argument of the sawtooth functions. This was introduced by 
Meyer [Me] and Dieter [Dil], and brought to a solid ground by Radenracher [R]: For a, b £ N, x, y £ R, the 

Dedekind-Rademacher sum is defined by 


s(a,b;x,y) := ^ 


k + y 

a —-- x 


k mod b 

Note that there is no loss in restricting x and y to 0 < x, y < 1. 


k + y 


The ideas of Apostol and Rademacher can also be combined: Takacs [T] introduced a shift in Apostol’s sum 

( 7 ): 


E * 


k + y 


n i k + y 

Jd+i a —;- x 


k mod b 

This was further generalized by Halbritter [H] and later by Hall, Wilson, and Zagier [HWZ] to the gener¬ 
alized Dedekind-Rademacher sum, defined for a , 6, c, m, n £ N, x,y,z £ R, by 
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On the other hand, we can start with the cotangent representation (4) of the Dedekind sum to arrive at 
generalizations. The equivalent of the Dedekind-Rademacher sum in terms of cotangents was first defined 
by Meyer and Sczech [MS], motivated by the appearance of the classical Dedekind sum in topology [HZ], 
The version we state here was introduced by Dieter [Di2]: For a, b, c £ N, x,y,z € R, the cotangent sum 
is defined by 


/ 7 . 1 / k -f- z \ / k z \ 

c(a, o, c; x, y, z) := - y c °t 7r I a - x ) c °t 7r I - D ) 

° k mod c ' ' ' ' 


Here the sum is taken over all k mod b for which the summand is not singular. Dieter remarked in [Di2] that 
the cotangent sums include as special cases various modified Dedekind sums introduced by Berndt [Be, BD]. 
Most of them are inspired by the transformation properties of the logarithm of the classical theta-function. 
We list them here using Berndt’s notation; throughout, a and b denote relatively prime positive integers. In 
the first sum, a, (3 G N, and a -1 is defined through a _1 a = 1 mod b. 



6-1 

5(o,6) := ^(-l) fe+1+ [ afe / 6 l 

k =1 


.,(0,6):=x;(-l) l “ k/ ‘ l ((f)) 

6—1 

s 4 (o,6) := 

k =1 

S5 (o,6):=^(-. 



Yet another generalization of (4) we mention here is due to Zagier [Za]. Generalizing the topological consid¬ 
erations in [HZ] to arbitrary dimensions, he arrived naturally at the following expression: Let ai,...,a^ 6 N 
be relatively prime to ao € N. Define the higher-dimensional Dedekind sum as 


s(ao; ai,..., aa) 


(_l)d/2 ^ nkai 

- - - — > cot-- 

a 0 ^ a 0 


7r ka d 

■ ■ ■ cot - 

a 0 


We note that this sum vanishes if d is odd, since the cotangent is an odd function. Berndt noticed in [Be] 
that a version of the higher-dimensional Dedekind sum was already introduced by Carlitz [C2] via sawtooth 
functions. 


Our Dedekind cotangent sums (Definition 1) combine the cotangent and higher-dimensional Dedekind sum. 
This could have been done by only introducing a shift of the variable in each cotangent of the higher¬ 
dimensional Dedekind sum. The reason for additionally introducing cotangent derivatives is twofold: first, 
they appear in lattice point enumeration formulas for polyhedra [DR]. We will make use of this fact in section 
6 about the computability of the Dedekind cotangent sums. Second, the cotangent derivatives are essentially 
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the discrete Fourier transforms of the Bernoulli functions. That is, our definition can be seen as the higher¬ 
dimensional ‘cotangent equivalent’ to Apostol’s Dedekind Bernoulli sum. In fact, the Dedekind cotangent 
sums include as special cases all generalized Dedekind sums mentioned in this section with exception of the 
Dedekind-Rademaclier sum and its generalization by Hall, Wilson, and Zagier. Technically, these sums could 
be treated in the same manner. What makes a subtle difference is the shift by a real number in the argument 
of the Bernoulli functions: to handle these, we would have to work with the discrete Fouries series of the 
shifted Bernoulli functions, which turn out to be much less practical than the ones without a shift. 

To be more precise, we state the discrete Fourier series of the Bernoulli functions next. In analogy to (3), 
we have the straightforward 


Lemma 6 For m >2, 


B 


m 




+ m 



p -1 

53 cot (m_1) 
k =1 


e 2*kn/p 


□ 


These discrete Fourier expansions can be used, for example, to rewrite the Dedekind Bernoulli sums in terms 
of the Dedekind cotangent sums: 


Corollary 7 If a,b, c € N are pairwise relatively prime and m, n > 2 are integers with the same parity then 


s m , n (a;b,c) d = £ (E) (A 

k mod a 


a 


— mn 2 m+n a m+n-l 
def (-l)( m -«)/ 2 

= mn --—-- 

O m-f-n 


(fE) + 


n m-\-n— 1 



f a 

b 

C ^ 

c | 

m + n — 2 

n— 1 

m — 1 


V o 

0 

o ) 


BmBn 

m-\-n— 1 


a 


We note that the parity assumption on m and n is no restriction, since the sums vanish if m + n is odd. It 
is worth mentioning that a close relative of these sums, namely, 


U fk 


fc=l 


EM; t - 


kb 


' 2 m+1 c 


a—l 

E 1 

k =1 


cot 


El cotC 
a 


( 7r kb 
V a 


appears naturally in the study of plane partition enumeration [Al]. 
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Proof of Corollary 7: By Lemma 6, 


n ,n(o;6,c)= B m (^\B n (— 

k mod a ' ' 


m+n 


a 

a— 1 


= mn — 
V 2a . 


to — 
V 2a . 


_ _T cot (m - 1} f-?-) cot (n_1) ) e 2 " fe(jh+ic)/a 

k mod a j,Z=l ' ' ' ' 

p _ Q ~^ 

-L^r?. V > V > . / /l / \ o~riUAh/ a 


E E 


E Ecot (m_1) (y) e 27riWa 

' ' fc mod a j=l ' ' 

v 7 v ' k mod a 1=1 v 7 v ' 


( 8 ) 


Now we use the fact that to + n is even, a, 6, c are pairwise relatively prime, and 

E 2-7T inkja f ® if a 117- 

“I 0 else. 

k mod a v 

Hence the first sum in (8) vanishes unless jb + Ic is divisible by a, that is, l = —jbc -1 mod a, where 
mod a. The second and third sum in (8) disappear completely, and we are left with 

. , , (-1 )(m+n)/2 «+ (nj\ ,„_u { 7 rjbc~ 1 \ B m B n 

’ v ' (2 a) m+n V a J \ a J a m+n - 1 

(_l)( m +n)/2 ^ } (n _ 1} / njb\ , B m B n 

\m+n a m+n-l h l^ COt \ a ) \ a 

j=l 


c x c = 1 


= mn- 


m+n—1 


3 Proof of the reciprocity law 

Proof of Theorem 2. Consider the function 

d 

f( z ) = cot<"*> 7T (a,jZ — Zj) . 

j =o 

We integrate / along the simple rectangular path 

7 = [x + zy, x — iy, x + 1 — iy, x + 1 + iy, x + iy\ , 

where x and y are chosen such that 7 does not pass through any pole of /, and all poles z p of / have imaginary 
part |Im( 2 p )| < y. By the periodicity of the cotangent, the contributions of the two vertical segments of 7 
cancel each other. By definition of the cotangent, 

lim cot(x ± iy) = =F i , 

y—*oo 

lim cot^ m \x ± iy) = 0 
y —*00 


and therefore also 
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If all nrij = 0, we obtain 


for m > 0. Hence if any of the nij > 0, 

[ f(z) dz = 0 . 

J 7 

f f(z) dz = i d+1 - {-i) d+1 = i d+1 (1 + (~l) d ) . 

J 7 

This can be summarized in 

1 ( i d if all rrij = 0 and d is even 

2i) 1 f{z)dZ = \ 0 otherwise, 

or, by means of the residue theorem, 

Tr^Res (f(z),z p ) = 


i d if all nij = 0 and d is even 
0 otherwise. 


(9) 


Here the sum ranges over all poles z p inside 7 . It remains to compute their residues. By assumption, / 
has only simple poles. We will compute the residue at z p = , k € Z, the other residues are completely 

equivalent. We use the Laurent expansion of the cotangent 


/ \ 1 / k + Zo \ 1 . 

cot 7r (aoz — Zo) = - z -+ analytic part , 

na 0 V a 0 J 

and, more generally, 

( ) . . (—l) m °mo! ( k + zo\ * mo+1 ' ) 

cot o) 7T ( a 0 z - z 0 ) = 7 - +1 -2-+ analytic part . 

(7ra 0 ) mo+i \ a 0 J 

The other cotangents are analytic at this pole: for j > 0, 


( 10 ) 


cot*"*'* 7 T ( a,jZ — Zj ) = E co ^ m 0 +h) x ( a 


L\ 

u> 0 3 


k + zo 
i 

a 0 


- Zi \ z — 


k + z 0 


a 0 


( 11 ) 


Hence 


o (tt \ k + zo\ (—l) m °mo! \ ' 77 

z lH-Md= m O 


d a- 

J _ Am-j+l 


V 


7T a 


k + z 0 
1 do 


Since 7 has horizontal width 1, we have do poles of the form fc +~° inside 7 , where k runs through a complete 
set of residues modulo qq. This gives, by definition of the Dedekind cotangent sum, 


k + z 0 \ 1 


E Hes ( f(z), z = — ^ ) = — (— l) m °rn 0 ! E 

k mod ao V °° J * h ,^ d 


* 1 , 

*iH-Hd= m 0 




d 0 

c | mo 

Zo 


dl 

mi + h 

Zl 


ad 

m d + Id 
Zd 


The other residues are computed in the same way, and give with (9) the statement. 


□ 
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4 Variations on a theme 


The conditions on the parameters appearing in Theorem 2 are not crucial; however, without them the 
theorem would not be as easy to state. Indeed, the conditions simply ensure that all poles of the function 
f used in the proof are simple. We will now state some special cases of Theorem 2 in which we drop some 
of the conditions. Strictly speaking, these ‘corollaries’ are really corollaries of the proof of Theorem 2. The 
first of those cases is the classical Dedekind reciprocity law, Theorem 1 [De]. Recall that 


s (a, b) 


1 

4 


( b 

a 

1 \ 

0 

0 

° 

V 0 

0 

0 ) 


Proof of Theorem 1. To modify our proof of Theorem 2, we have to consider the function 


f(z) = cot(7 raz) cot(Trbz) cot(7rz) . 


The residues are computed as above and yield the classical Dedekind sums; the only difference is an additional 
pole of order three at z = 0 (we may choose our integration path 7 such that 0 is inside 7). Its residue is 
easily computed as 

From here we can proceed as before. □ 


It should be mentioned that a proof of Dedekind’s reciprocity law along these lines is given in [RG, p. 21]. 
In fact, it was this proof that motivated the proof of Theorem 2. 


The second special case is an identity equivalent to the reciprocity law for the Dedekincl Bernoulli sums 
[Apl, Cl, Mi]. The respective statement for them can be obtained through Corollary 7. 


Corollary 8 Let ao,cti,a 2 € N be pairwise relatively prime, and mo,mi,m 2 G No not all zero, such that 
mo + mi + m 2 is even. Then 


l\ 1 2 
a/a 9 




l l 42 >0 
1 1 +^2 =m 0 


a 1 a 2 

mi + h m2 + I2 
0 0 


+ (-l) mi m!! ^ 


a n u a 


Iq 42 

1 0 + l 2 =rn l 


0 (X 2 
W-W. 


+(-ir 2 m 2 ! ]T 


l 0 4i>o 

i 0 + i l =m 2 


a o a i 

Wh\ 



a 0 a 2 

mo + 1 0 m 2 + 1 2 


0 


0 


0 


ao ai 

mo + lo m 1 + h 


0 


= ^(ao,ai,a 2 ,m 0 ,mi,m 2 ) 
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where 


(j)(ao,a 1 ,a2,m 0 ,m 1 ,m 2 ) := (-4 ) (m ° +mi+m2)/2 . 

1 '(- 1 )”*° v- f rn 0 \ ^ 2 fcl ^ 2 fc 

\ 2 fci — 1 — m\) k\k 2 




a 


mo+l 

'0 


E 

2fc 2 >m 2 -|-l 

2(fei+fc 2 —l) = m o+ m l+m 2 


: 2 2/ci — 1—mi 2/c 2 — 1—m 2 


(-i) r 


„mi+l 


2^0— 2fc 2 >m 2 +l 

2(/e 0 +fc 2 -l) = m 0 +m 1 +m 2 


TO 1 \B 2 koB 2 k 2 _ 2k 0 -l-m 0 2k 2 -l- 

2fc 2 -l-m 2 / fcofe ° 


(- 1 ) 


m 2 


a, 


,m 2 + l 


2fe n >m 0 + l, 2fc 1 >m 1 +l 


/ to 2 \ B 2ko B 2k 

\2fco — 1 — m 0 / fcofci 


1 2fc 0 -l-mo /1 2fci-l-mi 
a 0 


+ 


2(fc 0 +fc 1 -l) = m 0 +m 1 +m 2 

2m o + m 1+ m 2+2jBmo+mi+ma+2 / ( — 1 )( rn '°+ rn ' 1 ^ rra2 )/2,7t o!mi! q,™o+ m 1 + l 

to 0 + rni + m 2 + 2 V (m 0 + mi + l)!a™ 0+1 a™ 1+1 

(_l)(mo+m2-mi)/ 2rno , TO2!a mo+m2+l (_i)(mi+m2-mo)/2 mi!ffl2 l a ™i+™2+l 

(too + to 2 + l)!a™ 0+1 a™ 2+1 


(mi + TO 2 + l)!ar +1 ar +1 

Again the parity restriction on too + to 1 + m 2 is no constraint, since otherwise all the sums vanish. 

Proof. As in the last proof, we use the function 

f(z) = cot^ 1710 ^(iraoz) cot^^iraiz) cot ( ' m2 ^(a 2 nz) . 

By the pairwise-prime condition, the poles of / are all simple with exception of the pole at z = 0 (again we 
may choose our integration path such that 0 is contained inside). The residues at z = k/ao, for example, 
yield, as in the proof of Theorem 2, 


o 0 -l 


aX 2 


£ R* (/«,* = A)£ IM 

k =1 x u/ h,h>o 1 z 

l l+ l 2 =rn 0 


a 1 a 2 

TO 1 + h to 2 + l 2 

0 0 


Similar expressions are obtained for the other non-zero poles of /. To get the residue at z = 0, we use 
the following expansion of the cotangent, which follows directly from the definition ( 6 ) of the Bernoulli 
polynomials [BS]: 

I+^(- 1 ) fc 22 fci? 2 fc_ 2fe — 1 


cot z = 


k> 1 


( 2 * 0 ! 


( 12 ) 


Hence 


from which we obtain 


cot 


(rn) r — 


(- 1 ) 


m m\ 


^m -\-1 


___ ( 1 \kr\2k— 1 D 

(-1) £ 2k_ z 2k-l-r 


2fe>m+l 


k (2k — 1 — to)! 


Res (f(z),z = 0) =- 4>(a 0 ,ai,a 2 ,m 0 ,m 1 ,m 2 ) . 

7 r 

Since not all too,toi,to 2 are zero, the sum of the residues of / vanishes as in the proof of Theorem 2, and 
the statement follows. □ 

Next we prove the reciprocity law for the cotangent sums [Di2]. Recall that 

c (a,b,c;x,y,z) = c 


( c 

a 

h \ 

0 

0 

0 

V * 

X 

y / 
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It is easy to see that the reciprocity law of Meyer and Sczech [MS] is a special case of 


Corollary 9 (Dieter) Let a, b, c £ N be pairwise relatively prime, and define A,B,C by 


Abe + Bca + Cab = 1 . 


Let x,y,z € R, not all integers, and set 

x' = cy — bz, y' = az — cx, z' = bx — ay . 


Finally, let 


6{x) 


1 if x £ Z, 
0 else. 


Then 


c (a, b , c; x , y, z) + c ( b , c, a; y, z, x) + c (c, a, 6; 2 , x, y) = — 1- tS(z') cot*' 1 ! n (Acx 1 

ab 


(Be + Cb)y') 


— -r-b(x') cot^ 1 ! 7r (-Ba?/ — (Ca + Ac)z’) — —S(y') cot (1 ! 7r ( Cbz' — (Ab + Ba)x') . 
be ac 


Proof. As in the proof of Theorem 2, we use the function 

f{w) = cot 7 r (aw — x) cot ir(bw — y) cot ^ t{cw — z) . 

If / has only simple poles, the statement is a direct special case of Theorem 2. Otherwise, suppose we have 
a double pole w p , that is, there exist integers m and n such that (for example), 

to + x n + y 


To compute the residue of / at this pole, we use once again the Laurent series (10) and (11) of the cotangent 
to obtain 

Res (f(w),w = w p ) = —— cot (1 ! n ( cw p — z) . 

■nab 

We obtain similar residues for the simple poles as before, and the residue theorem yields an identity. That 
this identity is equivalent to Dieter’s can be easily seen by following the remarks in [Di2] just before Theorem 
2.3. □ 


The last special case is one in general ‘dimension’, Zagier’s higher-dimensional Dedekind sums [Za]. Recall 


that 


s (a 0 ;a lt ... ,a d ) 


( -l) d/ 2 



( a 0 

ai ■ ■ 

• a d \ 

c | 

0 

0 • • 

■ 0 



0 • ■ 

• o ) 


Corollary 10 (Zagier) If ao,..., Od £ N are pairwise relatively prime then 


where 


d 

y j5 (a n -a 0 ,...,c^,...,a d ) = l- h(a 0 ,..., a d ) , 

n —0 


h(a 0 , ...,a d ) 



E 

fc 0 >--. fc d >0 
fc()H-1 -kd = d/2 


B 2 k 0 ■ ■ ■ B 2 k d 2fc 0 

( 2 fc 0 )! • • • (2k d )\ a ° 


a 


2 k d 
d 
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This description of h(ao, ■ ■ ■, a d ) in terms of Bernoulli numbers was first used by Berndt [Be]; however, it is 
easily seen to be equivalent to the version given by Zagier. It is interesting to note that h(ao, ■ ■ ■, ad) can be 
expressed in terms of Hirzebruch L-functions [Za]. 

Proof. This time consider the function 


f{z) = cot7rao~ ■ ■ • cot 7 ia d z . 

As before, / has simple poles aside from the pole at z = 0. The residues are calculated as usual, for example, 


E Res (/(*)’*=£:) =^r E cot 


a o / 7ra 0 


irkai 

ao 


■ cot 


nka d (-l ) d / 2 


k =1 x 7 k =1 

The residue at 0 can be computed through rewriting (12) as 

\k o2fci 


ao 


s(a 0 ;ai,...,a d ) . 


( 1 \k rx2k TD 

:cotz = Y. ' - 


k>0 


(2 k)\ 


Hence 


m = - n V ■ 


o 0 ”-a d (7r^+ 1 , =0 - 0 


( 2 ^' 


and we obtain the residue 


Res(/(z),z = 0) 


^_l)d/2 2 d 

7 ra 0 • • • a d 


E 


£? 2 fc 0 ' ' ' I? 2 fe, 


-|-fc <i =d/2 


( 2 fc 0 )! • • • ( 2 fcd)! 0 


d 2fc 0 2/c d 


a A — 


(_i)rf / 2 


/i(a 0 , ...,a d ) . 


It remains to apply the residue theorem and (9). 


□ 


5 Petersson-Knopp identities 


Knopp applied in [Kno] Hecke operators to log 77 to arrive at Theorem 3. This identity was stated by 
Petersson in the 1970’s with additional congruence restrictions on a and b. For n prime, the Petersson- 
Knopp identity was already known to Dcdekind [De]. Theorem 3 was generalized by Parson and Rosen to 
Dedekind-Bernoulli sums [PR], by Apostol and Vu to their ‘sums of Dedekind type’ [AV], and, most broadly, 
by Zheng to what we will call sums of Dedekincl type with weight (mi, m 2 ) [Zh]. We will state and further 
generalize Zheng’s Petersson-Knopp identity after the following 


Definition 2 Let a, a \,..., a d € N. The sum 

S (a; Oi,..., a d ) := ^ /1 (—)••• f d (—) 

k mod a \ a / \ a J 

is said to be of Dedekind type with weight (mi,..., m d ) if for all j = 1,..., d, fj(x + 1) = fj{x) and 
for all a G N, 

fi 

k mod a 




a 


a mi fj{ax) . 


(13) 
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Note that the Bernoulli functions B m (x) satisfy (13) (with ‘weight’ — m+ 1), as do the functions cot^ m ^(7ra;) 
(with ‘weight ’m + 1). Zheng’s theorem is the ‘two-dimensional’ (d = 2) case of the following 


Theorem 11 Let n, a, a ±,..., a d € N. If 


S(a;ai,...,a d ) := ^ fi (—- 


k mod a 

is of Dedekind type with weight (mi,..., m<j) then 


ka-i 


, ,'ka d 
■ fd - 

’ a 


y ' t,~ mi - md 

b\n ri,...,r*d mod b 


E 


ai+na , ... , -ad + r d a 
b 


= n (Td^i-m 1 - m d [n) S (a; ai,..., a d ) 


Our proof is a relatively straightforward extension of Zheng’s proof for d = 2 [Zh]. We need the following 
two identities: 


Lemma 12 Let a, a \,..., a d G N. If 


s (a; T 


k mod a 

is of Dedekind type with weight (mi,..., m d ) t/ien for all j = 1,..., d, 




k mod a 


and 


S ( ab ; aid,..., a d 6) = b S (a; ai,..., a d ) 


(14) 

(15) 


Proof. (14): If (a, 6) = 1, the statement is vacuous. If a and b are not relatively prime, let a 1 = a/(a,b) and 
b' = d/(a, b). Then 


E u 

k mod a 



E E /, (•'• • k ~P) 

k mod a ' ' fc mod a' ' ' 

(a. &) E h ( X+ ^) = ^ b ) a ' mi = ( a > b) 1 ~ rrlj a m fj 

k mod a' ^ / 



(15): 


S' (ad; aid,..., add) 


E A 

k mod ab 




b E A 

k mod a 




b S (a; ai, ...,a d ) ■ 


□ 


Proof of Theorem 11. We will make use of two properties of the Mobius /i-function 


^(n) = 


1 

(-l) r 

0 


if n = 1, 

if n = pi • • - pm is square-free, 
otherwise, 
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namely, 


E b l ( d ) = 

d\n 


1 if n = 1, 

0 otherwise, 


ab ab/t 

E /( fc ) = E E ■ 

*=1_ t|6 fc=l 


These suffice to prove our statement: 

Y^ b -m 1 — -m d s(ab ; 


^ 6 _ roi _..._ ro<i ^ 5 (^afe; ■ 

6|n ri,...,rd mod 6 

= E^-- rorf E /d^ + ^j 

b|n mod 6 ' 

fc mod ab 

( = } E 6 _mi md E 

b\n k mod ab ' 

= E E (*,«— 

6|n k mod ab \ v / / 

= E /.fe) - 

bln clb mod ai> 

tt- M — 


nai nad 

— -h rya,..., —— + r d a 

b b 


E (knar , , (kna d , kr d 

ix 1 + — )'"J d { + — 


kna\ 
ab(a, b ) 


ab 2 b J 
{k, b) 1 ~ md b md f d 


kna d 
ab(a, b) 


k mod ab 
(k,b)=c 


knar 

ab(a , 6) 

knai \ 


/ kna d 
\ab(a, b) 

kna d \ 


= EE c 

6|n c|fe 

( =EE< 

b|n c|6 


,d-mi- rrid 


i—m 1 - 


EE 

m d E n(t) S ( ■ 

b\n c\b 

t\b/c 

E c d mi - 

S(ae;^, 
V c 

c£e|n 

E^-- 

ce|n 

- md S (ae; 

V c 

> ^ c d- mi -. 
ce=n 

' _m<i S' (ae; aie,..., 

) 

Tl (7d—\—m\ — 

— main) S (a;ai,... 


(fc,b/c) = l 

‘E*o E /.(^)-a(^ 

£|b/c /c mod ab/ct 

\T^ /,N O / ^1 nC I d \ 


The fact that the Dedekind cotangent sums 

( ao a\ 
mo mi 
0 0 


Z)= e ncot(-)^ 

0 / k mod ag J—1 


have weight (mi + 1,..., rrid + 1) immediately yields the Petersson-Knopp-like Theorem 4. 














Dedekind cotangent sums 


15 


A particularly simple form of Theorem 4 is achieved for Zagier’s higher-dimensional Dedekind sums (the 
case mo = • • • = rrid = 0): with 

n cr_i(n) = E = E d = CT ( n ) > 

d|n d|n 


we obtain 


Corollary 13 For n, ao, ■ .., ad € N, 

^ s(a 0 6;-^ai + ria 0 ,..., + r d a 0 ) = cr(n) s(a 0 ; ai,..., a d ) . 

6|n ri,...,r d mod 6 


□ 


6 Proof of the polynomial-time computability 


In this last section, we prove that the Dcdekind cotangent sums are computable in polynomial time (Theorem 
5). Our proof is similar to the last section of [BR], where the computability of Zagier’s higher-dimensional 
Dedekind sums was shown. Again we will ‘merge’ two theorems in combinatorial geometry, due to Barvi¬ 
nok [Ba] and Diaz-Robins [DR]. The latter allows us to express generating functions for the lattice point 
count in cones in terms of cotangents; Barvinok’s theorem tells us that the respective rational functions are 
computable. The usual generating function for lattice point enumeration for a d-dimensional cone /C, 

F K (< l ) = F K (q u ...,q d ):= £ q?' • • • q?‘ , 

m=(mi,...,md)(zJCr\Z d 


assumes in [DR] an exponential variable: 


Theorem 14 (Diaz-Robins) Suppose the cone K. cM. d is generated by the positive real span of the integer 
vectors Vi ,...,v d € Z d such that the d x d-matrix M = (), whose column vectors are vi,...,v<j, is 
lower-triangidar. Let pk '■= an • • • akk (k = 1,..,, d) and G := (Z/piZ) x • • • x (Z/p^Z). Then 


E 

iG/cnz d 


• 27 r(m.s) _ 


2 d |G| 


En 

reGfe=l 


1 + coth — (s + ir, Vfc) 
Pk 


Here ( , ) denotes the usual scalar product in . 


We note that the assumption on M being lower-triangular is not crucial for practical purposes: by Hermite 
normal form, any cone generated by integer vectors is unimodular equivalent to a cone described by a 
lower-triangular matrix. 


Theorem 15 (Barvinok) For fixed dimension d, the rational function F/c( q) is polynomial-time com¬ 
putable in the input size of 1C. 


These two powerful theorems combined allow us to compute the Dedekind cotangent sum in polynomial 
time: 
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Proof of Theorem 5. It suffices to prove the polynomial-time computability of 


d 

ri cot ^ 

k mod a j —1 



We will do this inductively by constructing a cone whose generating function has the above Dedekind 
cotangent sum + ‘lower-dimensional’ sums, that is, Dedekind cotangent sums with a lower number of factors. 
More precisely, we will construct a cone whose generating function will be 



By induction and Barvinok’s Theorem 15, this will prove the computability of the Dedekind cotangent sum. 
Let K, be the positive real span of 


( a \ 

0 


0 

\hj 


( 0 \ 

a 

0 

0 

\bij 


( o \ 


0 
a 

V b d~l ) 


( 0 ^ 

0 

V a J 


The integers b ±,..., bd-i are to be chosen later. We will repeatedly use the fact that 


coth 7r 

k mod a 



acothiraz . 


(18) 
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To apply Theorem 14 to /C, note that in our case pk = a k . Hence 


E « 


—2-7r{m,s) 


1 


2d a d(d+l)/2 


a _ 

e n ( i+c ° th j( s +* r ’ vfc 


mod aO k— 1 


1 

2d a d(d+l)/2 


Y (l + c °th ^ ((si + iri)a + (s d + ird)&i)) 


mod aO 

d=i,....d) 


(l + coth ((s 2 + ir 2 )a + (s d + ir d )fo)) ■ • • (l + coth -^( s d + ir d )aj 
— E + C0th7T + {s d + ir d)^yj 


2d a d(d+l)/2 


mod aO 

0=3,..., d) 


(l + coth ((s 3 + ir 3 )a + (s d + ir d )b 3 )) • • • (l + coth -^y (s d + 


^ | 1 + coth - 

r 2 mod a 2 


+ — + (s d + ir d )- 


( 18 ) 


1 


2d a d(d+l)/2-3 


E ^l + coth7r^Si + (s d + *r d )h^ 


mod aJ 
0=4,...,d) 


1 + coth 7T ^s 2 + (Srf + ir d )y^j ^ (l + coth ((S 4 + *r 4 )a + (s d + 

• ^l + coth-^ T (s d + «r d )^ ^1 + coth?r + (s d + *r d )h 


7*3 mod a 3 


(18) 


2 d a d 


Y ^1 + C0th7T + (s d +ir d)~ S j S j ' ' 


mod a d 


1 + coth 7r ^s d _i + (s d + ir d )E^^ (l + coth-^ T (s d + ir d )J 

,EE( 1 + coth 7i ^si + ( s d + i(na + k )) — ^ ' 


1 + coth7r ^s d -i + (s d + i(na + k)) (l + coth ^Ef ( Sd + *( na + &))) 

^ ^1 +coth7r ^si + (s d + ik )~^ ''' 


1 + coth 7r ( s d -i + (sd + ik) ¥))£( 1 + coth 


zn 


Sd , ik 


H—:r^r + 


(18) 1 


fc=i 


= 2 ^ E ( 1 + coth7r ( si + (*<*+ i/c )- 




1 + coth7r ( Sd-i + (sd + ] ] [ 1 + coth7r ( — + 


Sd ik 
a a 


If we now choose bj = a d l aj (j = 1,..., d — 1), where a d L a d = 1 mod a, and 

Sj = - iz d bj (j = 1,..., d - 1), Sd = ia^d , 
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this generating function becomes 



^1 + COt 7T 



□ 


7 Closing remarks 

There remain many open questions. First, for all Dedekind sums there are closed formulas for special cases, 
for example, s(l,a). Many of such formulas can be found in [BY]. The Dedekind cotangent sum and its 
variations attain such closed formulas for special variables; it is not clear how far these cases lead. Second, 
it seems desirable to have a parallel theory for the respective Dedekind sum in which the cotangents get 
replaced by Bernoulli functions. This will most certainly require different methods than the ones used in 
this paper. Finally, we can apply our very general Petersson-Knopp-like Theorem 11 to various other sums, 
for example, to such 'Dedekind Bernoulli sums’. 

Acknowledgements. I am grateful to Sinai Robins for good discussions and Thomas Zaslavsky for many 
helpful comments on a previous version of this paper. 
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